The Peierls substitution in an engineered lattice potential 
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Artificial gauge fields open new possibilities to realize quantum many-body systems with ultracold 
atoms, by engineering Hamiltonians usually associated with electronic systems. In the presence of 
a periodic potential, artificial gauge fields may bring ultracold atoms closer to the quantum Hall 
regime. Here, we describe a one-dimensional lattice derived purely from effective Zeeman-shifts 
resulting from a combination of Raman coupling and radiofrequency magnetic fields. In this lattice, 
the tunneling matrix element is generally complex. We control both the amplitude and the phase 
of this tunneling parameter, experimentally realizing the Peierls substitution for ultracold neutral 
atoms. 



Ultracold atoms subjected to artificial gauge fields can 
realize phenomena usually in the domain of electronic 
systems. Prime examples include the quantum Hall effect 
(for abelian gauge fields), and topological insulators (for 
non-abelian gauge fields) PQ . Many of these phenomena 
are predicted to occur at extremely low temperatures, 
and adding a lattice potential to an ultracold system can 
increase the energy scales at which strongly correlated 
states are expected to emerge [H [5] . Current techniques 
for generating periodic potentials in ultracold atom sys- 
tems use optical standing waves created with suitably 
polarized counterpropagating lasers [3]. In contrast, we 
describe a one-dimensional (ID) "Zeeman lattice" for ul- 
tracold atoms created with a combination of radiofre- 
quency (rf) and optical-Raman coupling fields, without 
any optical standing waves. In this lattice, atoms acquire 
a quantum mechanical phase as they hop from site to site, 
explicitly realizing the Peierls transformation [5] in the 
laboratory frame. Our approach extends existing Raman 
dressing schemes [5] by simultaneously generating an ar- 
tificial gauge field and an effective lattice potential. 

Optical lattices generally result from the electric dipole 
interaction between an atom and the electric field of 
an optical standing wave, yielding a potential Vdi p (r)cx 
a(\)I(r), where a(X) is the atomic polarizability at laser 
wavelength A, and I(r) is the spatial intensity distribu- 
tion pQ. In such lattices, the natural units of momen- 
tum and energy are given by the single photon recoil 
momentum hkr, = 2nh/X and its corresponding energy 
Er 1 = h 2 k 2 ^/2m, where m is the atomic mass. 

Quantum particles with charge q in a ID periodic po- 
tential (here along e^) acquire a phase (f>j = (q/ti) J 3+1 A- 
e x dx upon tunneling from site j to j + 1 in the presence 
of a vector potential A. For sufficiently strong potentials, 
this system is described by the tight-binding Hamiltonian 

H = -^[iexp(i^-)a] +1 % +h.c], (1) 

where describes the creation of a particle at site j, 
and texp(i(j)j) is the complex matrix element for tun- 



neling between neighboring sites. Using the phases <j>j 
to represent the effect of A is known as the Peierls 
substitution [5], and for uniform phase <fi the energy is 
E(k x ) = — 2t cos{ / Kk x /kr l — <fr), where k x is the particle's 
crystal momentum. 

We realize the Peierls substitution for ultracold atoms 
by synthesizing a ID effective Zeeman lattice that allows 
independent control of both t and <p. Previous experi- 
ments (a) controlled the amplitude and sign of t in driven 
optical lattices [7j, or in addition (b) controlled cf> by 
means of rotating optical lattices [3] or Raman-assisted 
tunneling in an optical superlattice [5] . Our effective Zee- 
man lattice technique provides both a periodic potential 
and an artificial vector potential in the laboratory frame. 

The Zeeman lattice arises from a combination of rf and 
Raman fields that simultaneously couple the spin states 
{\mF)}m F =o,±i of 87 Rb's F= 1 ground level, which are 
split by hwz (Fig. [lp.,b). In the frame rotating at the rf 
frequency Au) and under the rotating wave approxima- 
tion, the combined rf-Raman coupling contributes a term 

H l{+R (x) = ^(x)-F + H Q (2) 

to the overall Hamiltonian, where F = (F x ,F y ,F z ) 
is the F — 1 angular momentum operator; Q = 
[n T f+fl^cos(2k[ J x)— riR,sin(2fcia;),v / 25) /V2, in which 
f2 r f and S1r are the rf and Raman coupling strengths, 
and 5 = Auj — ujz is the detuning from Raman resonance; 
and HQ=—e(h 2 I—Fz)/h describes the quadratic Zeeman 
shift. Equation |2]) is the Zeeman Hamiltonian for an 
effective field B G ff (x)=htt(x) / gF [J-Bi where fiB is Bohr's 
magneton and qf is the Lande g-factor. This spatially 
varying effective Zeeman shift produces a ID lattice po- 
tential (Fig. [lp). As atoms tunnel from site to site, B c g 
rotates by 2-7T about e z (Fig. fljl) and the atoms acquire 
a geometrical Berry's phase jTOj proportional to the en- 
closed solid angle Ob- The tunneling parameters t and 
</>, obtained from E(k x ), are non-trivial functions of ft. 

When O r f SIr., <5, e the effective Zeeman shift re- 
duces to h\Cl\ h[Q r f + J7r cos(2kLx)]/^/2, and when 
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FIG. 1. Effective Zeeman lattice, a-b. A uniform magnetic 
field B e y Zeeman-splits the levels in the F — 1 ground state 
manifold of S7 Rb by cjz, and provides a quadratic Zeeman 
shift e. In conjunction with an rf magnetic field B T fe x with 
frequency Acj, a pair of orthogonally-polarized counterpropa- 
gating Raman beams with frequencies (uj,uj+Auj) illuminates 
the atomic sample. The rf and Raman fields have coupling 
strengths Q r f and Or. c. The spatially varying eigenvalues of 
Hii+n(x) (red, blue and green curves) give rise to our A/2 
effective Zeeman lattice; as plotted Ml I { = lEi, hSln = 10EL, 
and H8 = 2El. d. Spatial precession of B c ff (x) (dark arrow) 
and the solid angle 0b it subtends when an atom tunnels 
to the nearest neighboring site (points 1 to 5 in c). This 
geometrical Berry's phase gives the Peierls phase <j>. 



rip- 3> fi r f,<5 we obtain the analogous result H\il\ f=s 
^[fijj, + O r f cos(2/c/,a;)]/v / 2. In both of these limits, the 
larger of the two fields defines a natural quantizing axis 
about which the smaller field spatially modulates |fi|. 
For rifi^firf, this quantizing axis is spatially rotating. 

We experimentally characterize the lattice in three 
ways: (i) we measure the effective mass m* — 
h 2 [d 2 E(k x )/dk x ]~ 1 , which in the tight-binding regime is 
inversely proportional to f; (ii) we quantify the Peierls 
phase (j> and test its robustness against fluctuations in 
firf ; and (iii) we investigate the diffraction of BECs from 
our effective Zeeman lattice. In each case, we start 
with 87 Rb BECs in the \F = \m F = -1) state in a 
crossed optical dipole trap with frequencies (f x , f y , f z ) = 
(13,45,90) Hz [TT]. In the presence of a uniform bias 
field Boe y , we apply an rf magnetic field with fre- 
quency Aw/27r= gru^Bo =3.25 MHz and prepare the 
BEC in the lowest energy rf-dressed state [12]. Two 
A = 790.33 nm Raman laser beams, counter-propagating 
along and differing in frequency by Aw, couple the 
BECs internal degrees of freedom with strength 
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FIG. 2. Effective mass. a. Comparison of the oscillations of 
a BEC in the \mp = —1) state to those in an rf-Raman 
dressed BEC [M2 R = 12.4(9) E L and Kl xi = 2.04(6) E L ]. The 
curves are fits to a sinusoid from which we obtain f x = 
14.0(1) Hz and /* = 5.3(1) Hz, thus m*/m = 7.0(3) and 
t — 0.015(1)El- b. Measurements of m* /m as a function of 
S7r and fl T {. The curves depict the expected m*/m ratio. 

(Fig. [jj,,b). The combination of rf and Raman coupling 
creates a ID lattice potential along e x , the direction of 
momentum exchange defined by the Raman beams. 

We obtain the atoms' effective mass m* by inducing 
dipole oscillations [T3HT5] along and measuring shifts 
in the oscillation frequency as a function of the coupling 
strengths and f2 r f. The atoms slosh in the lattice 
for a variable time r, after which we remove all coupling 
and confining potentials (thus projecting the final spin- 
momentum superposition into bare atomic states) and 
absorption image after a 28.2 ms time-of-flight (TOF). 
Figure [2^, shows bare and dressed condensates oscillat- 
ing at frequencies f x and /*, respectively; Fig.(2}} shows 
that the effective to bare mass ratio m* /m = (f x /f*) 2 
as a function of JIr and O r f is in good agreement with 
calculations (curves [Hj). These data provide the tun- 
neling matrix element amplitude t/E^ = (m/m*)/ir 2 in 
the tight-binding regime |16j . 

An important characteristic of our effective Zeeman 
lattice is the presence of a tunable Peierls hopping phase 
</>, which can be revealed through its effects on E(k x ) and 
is experimentally controlled by adjusting fl z . We quan- 
tify (f> both by adiabatically modifying the bandstructure 
( "adiabatic method" ) and by inducing oscillations in mo- 
mentum space ( "sudden method" , similar to above). Fur- 
thermore, we test its robustness to variations in J7 r f with 
the latter method. 

In the adiabatic method, wc load a BEC at k x = 
and adiabatically change Q z , such that the BEC always 
sits at the minimum of E(k x ) located at A: m i n . The time 
scale for adiabaticity is set by the modified trapping fre- 
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FIG. 3. Peierls transformation, a. Peierls phase <fi measured 
using adiabatic (crosses) and sudden (circles) changes of f2 z . 
Vertical lines denote the first Brillouin zone. b. Sloshing 
amplitude after suddenly changing Q. z . We observed strong 
damping of oscillations in the region shaded in gray. c. Tun- 
neling amplitude t measured from oscillation frequency. The 
rf coupling was modulated as a function of Q, z to test the 
robustness of the Peierls phase <j>. The Raman coupling was 
held at ftf^R, = 10.0(8)i?L. The dashed curves correspond to 
the expected behavior calculated from i/ r f + a, and the pink 
bands arise from the experimental uncertainty in Or. 

quency /* along the direction of the Raman beams. Once 
Q z reaches its final value, we remove the trapping poten- 
tial and subsequently deload all atoms into the |mi? = + l) 
spin state while mapping the occupied crystal momentum 
k x to free-particle momentum |I4) . We image this distri- 
bution after a 13.1 ms TOF, revealing k m i n . The Peierls 
phase, shown as crosses in Fig. [3]a,, is <fi/ir = fc m in/&L- 

In the sudden method, we test the robustness of the 
Peierls phase 4> by first adiabatically loading to <f> = ±ir 
(the condensate sits at the edge of the Brillouin zone) 
and then suddenly changing both Q, z and f2 r f [17] to 
new values (changing both ip and t). This results in 
momentum space oscillations centered at fc m i n . After 
a time r we release the BEC, and measure as above. 
We fit the crystal momentum dynamics with k x (r) — 
k min + Ak x cos(27tt/* + 7), where Ak x is the amplitude, 
and 7 is an overall phase-shift whose average value is 
0.9(l)7r for these measurements. Figure[3^ (circles) shows 
the measured Peierls tunneling phase as a function of Q z . 

Measurements from the adiabatic and sudden methods 
are in good agreement with each other and their expected 
values (Fig. [3^,, dashed curves), highlighting the precise 
experimental control offered by our rf-Raman induced 
effective Zeeman lattice. This agreement also demon- 
strates the robustness of our engineered Hamiltonian to 
deliberate variations in Q r f of up to 0.25El. 



The sloshing amplitude \Ak x \ is displayed in Fig. [3(3. 
For large initial \Ak x \ (shaded region) the oscillations are 
strongly damped, which we attribute to energetic insta- 
bilities resulting in depletion of BEC atoms [IS] ■ This is 
evident from the departure of the oscillation amplitude 
from the expected value; the experimental range of the 
dynamical instability, h\fl z \ < I.IEl (corresponding to 
£^k x > 0.6fc^), is shaded in gray and is in agreement 
with previous observations in conventional optical lat- 
tices [IS]. Figure [3fc displays the tunneling amplitude t, 
obtained from /*. For comparison, a sinusoidal lattice 
would require a depth Vq ~ 8El to give similar parame- 
ters. 

Having discussed the behavior of atoms in the lat- 
tice's lowest band, we now explore the full lattice by 
suddenly turning it on, diabatically projecting a ground 
state BEC into higher bands. At the beginning of such a 
pulse(r pu i so <C irh/y/sEL, where s = Vq/El), an ordinary 
periodic potential, would spatially modulate the BEC's 
phase [H|; our effective Zeeman lattice induces such a 
modulation but in a spin-dependent manner. We focus 
on the Or 3> fi r f and £^r <C O r f tight-binding regimes 
and investigate the spin and spatial structure of our lat- 
tice. Our data extends well beyond the short-time phase 
modulation regime. 

In the absence of either Raman or rf coupling, there is 
no lattice. As indicated in Fig. [4^,, we use two different 
methods to introduce our lattice on an initial spatially 
uniform state: (i) starting with an rf-dressed state (with 
k x =0), we suddenly (< 1 /us) turn on the Raman beams; 
or (ii) starting with a Raman-dressed state [12] (a super- 
position of \mp = 0, k x — 0) and \mp =±1, k x =p 2ki)), we 
suddenly turn on the rf-ficld. 

After holding the lattice on for a time r pu i se , we sud- 
denly turn off the rf and Raman fields, together with the 
confining potential. The atoms are projected onto the 
bare spin-momentum basis and separate in TOF in the 
presence of a magnetic field gradient (along e^), allowing 
us to resolve their spin and momentum components. 

We observe detectable population in states with mo- 
menta up to \k x \ < AkL (Fig. |4}d). We perform such 
experiments for f2R,/f2 r f sa 3 and 5. We minimize 
the effects of interactions by working with small BECs 
(?« 9 x 10 4 atoms). Figures [4j:,d show the fraction of 
atoms in each diffracted order evolving with time. We 
observe multiple revivals of the initial spin-momentum 
state and find symmetry in the population dynamics of 
spin-momentum states with opposite momentum and op- 
posite spin. The curves represent fits to the popula- 
tions in all spin- momentum components. The parame- 
ters from the fits are all within 10% for our calibrated 
values, demonstrates that the spin-momentum dynamics 
are well described by the unitary evolution of the initial 
states under -ff r f+R [14] - 

Based on this technique for controlling the Peierls 
phase and inspired by recent proposals for creating flux 
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FIG. 4. BEC diffraction from the effective Zeeman lattice, a. Starting with an rf-dressed (Raman-dressed) state, we suddenly 
turn-on the Raman (rf) field for a variable time r pu i sc . b. Using TOF absorption images of the projected spin-momentum 
distributions, we count the number of atoms in each diffracted order and determine its fractional population. Panels c,d depict 
time evolution of these fractions. The curves are fits to the data, calculated from H t { + r. The fit parameters are: (c.) rf-dressed 
h(Qrf, Or, fi z )=(3. 57, 11.49, -0.04)£ L and (d.) Raman-dressed h(Q r! , fi R , fi z )=(3.06, 15.14, 0.0S)E L . 



lattices [TJ [20] > we now describe how this method might 
be extended to create a lattice with zero net flux that is 
topologically equivalent to the Hofstadter model with flux 
density n$ = 1/3 per plaquette. Because the hopping 
phase is only denned modulo 2-7T (thus n$ is only defined 
modulo 1), a uniform magnetic field with n$ = 1/2 is 
equivalent to a staggered field with n$ = ±1/2. In the 
same spirit, a magnetic field staggered along e y with flux 
density (. . . , 1/3, 1/3, —2/3, . . .) has zero net flux yet is 
equivalent to a uniform field with n$ = 1/3. These fields 
could be generated by the Peierls phases 4> y {jx,jy) = 
and 4> x (j x , j y ) = — (2tt/3) mod (j y , 3) . Reminiscent of the 
flux rectification mechanism proposed in Ref. |21| . this 
configuration can be created in our system by adding two 
standing waves along e y (normal to the Raman lasers): 
a state-independent lattice with period a localizing the 
atoms to specific lattice sites (e.g. from a retro-reflected 
532 nm laser), and a state-dependent "vector" lattice 
with a period 3a/2 sinusoidally modulating £! 2 as a func- 
tion of y (e.g., from an additional 790 nm laser, nearly 
counter-propagating). Figure^ shows that with a suit- 
able relative phase between the standing waves, the fluxes 
along e y repeat with the pattern (—2/3,4/3, —2/3), giv- 
ing the desired flux per plaquette. To verify this heuristic 
interpretation, we numerically solve the 2D bandstruc- 
ture (for exact parameters see [14]). and as shown in 
Fig. [5)3, we confirm that for a wide range of parameters, 
the three lowest bands are described by the same Chern 
numbers [22] as are those of the n$ = 1/3 Hofstadter 



model: (1,-2,1). 

We realized a ID lattice potential for ultracold atoms 
using only rf and Raman transitions, in which the tunnel- 
ing matrix element is in general complex. This work con- 
stitutes a first step towards realizing flux lattices [20] , i n 
which the physics of charged particles in strong magnetic 
fields can be simulated. The tunability of the Peierls 
phase achieved with our rf- Raman lattice would allow the 
observation of nonlinear effects of ultracold atoms in ID 
periodic potentials, such as atomic density modulations 
with periodicity larger than the lattice spacing 23J. 
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tier Center. K.J.-G. thanks CONACYT; L.J.L. thanks 
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Supplementary Material 
INDUCING OSCILLATIONS 

Starting with the rf-dressed state described in the main 
text, we ramped the Raman beams from to Or > 4 
El /fi in 70 ms. Then we ramped the rf coupling strength 
to an adjustable final value f2 r f in 2 ms, such that an 
effective Zeeman lattice was created. We induce sloshing 
by applying a synthetic electric field [15] . achieved by 
ramping the z-component of the Zeeman field to Ml z ~ 
2El and then jumping it back to in 2 ms. 



MOMENTUM REPRESENTATION AND 
DIAGONALIZATION OF H rf+R {x) 

In order to compute the properties of our lattice poten- 
tial we numerically diagonalized the combined rf-Raman 
Hamiltonian -f/rf+R. given the experimental parameters 
f2 r f , JIr, Q z . For simplicity we work in the momentum 
space representation of ifrf+R. This representation offers 
an alternative understanding of the effective Zeeman lat- 
tice structure that arises from the combination of rf and 
Raman coupling fields. 

While Raman transitions couple together states with 
mp differing by ±1 and momentum differing by ±2hkL, 
rf-coupling processes only change mp by ±1, leaving the 
momentum unchanged. The combination of both Raman 
and rf fields sets appropriate conditions for secondary 
rf and Raman processes to populate states with higher 
momentum. The available states under rf-Raman cou- 
pling constitute a set of spin-momentum states = 
{\rriF, ti{k x + 2nfci))} where n G Z; mp = 0, ±1. As ex- 
pected, this basis is that of a lattice. 

In the basis {|\&n)}, -f/rf+R, is a Hermitian block ma- 
trix of size i{2N J \-\), when n is restricted to n < N. For 
our parameters, dimensions larger than 3(2A^ + 1) = 81 
provided indistinguishable results. The dimension of this 
basis is appropriate for our calculations since we observed 
the population of states with up to |n| =4. For clar- 
ity, we construct the Hamiltonian by arranging the spin- 
momentum states in the {^n} basis with increasing mo- 
mentum index n = — N,— N+l, . . . ,N. Along the princi- 
pal diagonal, we have 3x3 blocks 

A K (n)=h 2 {k x +2nk L ) 2 I+ ^F x ~ [n z F z + ^(h 2 I-F 2 )]; 

these terms correspond to kinetic energy, rf coupling of 
spin states with equal momentum, and the real Zeeman 
interaction, respectively; {F x , F y ,F z } are the matrix rep- 
resentations of the F = 1 angular momentum operators, 
and I is the 3x3 identity matrix. Above and below the 
A fcx (n) blocks, we have 3x3 blocks B = y/2£l- R (F x —iF y )/4: 
describing the Raman coupling of spin-momentum states 
differing in momentum by Ak x = ±2k.L. 
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The diagonalization of i? r f + R,(n) as a function of k x 
gives the bandstructure of the combined rf-Raman lattice 
potential, E(k x )—E(k x , f2 r f, Or, fl z , <fi). 



by computing the shifts the lowest band [with minimum 
at fc m j n = (0/7r)fci] as a function of the experimental pa- 
rameters. 



EFFECTIVE ZEEMAN LATTICE PROPERTIES 

We extract the Zeeman lattice properties by fully char- 
acterizing the lowest energy band which in the tight bind- 
ing regime is of the form E(k x )—— 2t cos(irk x /kL—(fi). The 
tunneling amplitude is given by t — AE/A, where AE is 
the width of the lowest band. The lattice depth was cal- 
culated as the depth that a lattice potential would have 
in order to produce the calculated width AE. The effec- 
tive mass is defined as m* — h 2 [d 2 E(k x )/dk x ]~ 1 , where 
the derivative is evaluated at the point of interest (e.g. 
for our Q z = measurements of effective mass displayed 
in Fig. 2, this was at k x — 0). In the tight-binding 
regime, t and m*/m are inversely proportional to each 
other n 2 t/EL =m/m*. We obtained the Peierls phase <j> 



DELOADING 

We deloaded to a single bare spin state by rapidly 
ramping to Q z — 0, taking ft£! R — > to zero in 500 fis 
while increasing M7 r f — > 3El. We then ramped the 
M7 Z — 140El transferring all atoms into |tof = +1). 

2D LATTICE FOR 1/3 FLUX 
HOFSTADTER MODEL 

We numerically solve the bandstructure of the 2D lat- 
tice arising from the combination of a period a — 266 nm 
lattice with depth Vb = 22.5El, and a ID "Zeeman 
lattice" with O rf = IE L , R = 10E L , e = 0A4E L . 
Here El denotes the Raman recoil, not the short-period 
lattice recoil. 



